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The ground state of the S — ^ Heisenberg antiferromagnet on the pyrochlore lattice is theoret- 
ically investigated. Starting from the limit of isolated tetrahedra, I include interactions between 
the tetrahedra and obtain an effective model for the spin-singlet ground state multiplet by third- 
order perturbation. I determine its ground state using the mean-field approximation and found 
a dimerized state with a four-sublattice structure, which agrees with the proposal by Harris et 
al. I also discuss chirality correlations and spin correlations for this state. 
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Antiferromagnetic spin systems on geometrically frus- 
trated lattices show many_unusual behaviors in magnetic 
and thermal properties .ErH^ Their origin may essentially 
be attributed to the point that their classical ground 
state has a thermodynamic degeneracy with a poten- 
tial instability towards various phases. Villain pointed 
out the possibility that thermal fluctuations drive syrp.- 
metry breaking to the phase with the largest entropy.™ 
This mechanism of order by disorder is also expected for 
quantum fluctuations. In this paper, I study the ground 
state of the S — ^ Heisenberg antiferromagnet on the 
pyrochlore lattice and examine the possibility of a new 
quantum phase. 

The pyrochlore lattice is a three-dimensional (3D) net- 
work of corner-sharing tetrahedra (see Fig. 1), realized 
as a sublattice of many compounds including pyrochlorcs 
and spinels. Since all the elementary closed paths of spin 
exchange processes are triangular, antiferromagnets on 
the pyrochlore lattice are subject to strong frustration 
due to this lattice geometry. Another important charac- 
ter is the loose packing of the lattice, which is in contrast 
with the closed packing of the fee lattice, another exam- 
ple of a geometrically frustrated lattice in 3D. 

The combination of geometrical frustration and loose 
packing of the pyrochlore lattice leads to thermody- 
namically degenerate ground states of classical antifer- 
romagnets on it, which was first pointed out by An- 
derson for the Ising modelD^ Numerical simulations for 
the classical antiferromagnets indicate the absence of a 
finite-temperature magnetic phase transitionJiP suggest- 
ing that thermal fluctuations cannot break symmetry. 

Quantum fluctuations are another jphannel of lifting 
the classical ground-state degeneracy.EP This possibility 
for thfi pyrochlore lattice was first examined by Harris 
et al!& They used a field-theoretical approach and de- 
generate perturbation for the spin-singlet ground-state 
multiplet and showed a sublattice ordering of dimeriza- 
tion. This was followed by Canals and Lacroix, who 
also started from decoupled tetrahedra and calculated 
the spin-spin correlation function at finite temperatures, 
using perturbation in the tetrahedron-tetrahedron inter- 



actions up to the second order.© The result shows that 
the spin correlation length remains very short at the low- 
est temperature, and they concluded that the ground 
state is a quantum spin liquid with a finite spin gap, 
which is consistent with the claim of Harris et al. In 
the pyrochlore lattice, it is believed that there exist a 
thermodynamic number of singlet states below the spin 
gap. In this paper, I study the possibility of symmetry 
breaking in the singlet sector driven by quantum fluc- 
tuations. To this end, I follow a similar approach as 
those of the previous two works and start from a set of 
weakly interacting tetrahedra. I construct a full effective 
model for the spin-singlet ground-state multiplet within 
the third-order degenerate perturbation. I show that its 
mean-field solution is a partially dimerized state with a 
four-sublattice structure, and discuss spin and chirality 
correlations. 

The pyrochlore lattice consists of two types of tetra- 
hedra with the same number, say type A and B, and 
the network of these tetrahedra is bipartite as shown in 
Fig. 1. I study the Hamiltonian of the S — ^ Heisenberg 
antiferromagnet on this lattice, 




Fig. 1. Pyrochlore lattice and a unit of type-A tetrahedron. 
Tetrahedra of type A (B) are shown by black (white) bonds. The 
lattice constant of the cubic unit cell is set to units of length. 
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Fig. 2. Eigenenergies of the Heisenberg antiferromagnet on (a) a 
tetrahedron pair and (b) a tetrahedron triple. 

if = J ^ Sr • Sr' + J' X! • Sr-, (J, J' > 0)(1) 
(r,r')eA (r,r')eB 

where the first and second summations are taken over 
bonds on type-A and B tetrahedra, respectively. To 
study the case of homogeneous couphng J'— J, I start 
from the hmit of J'=0 and include J'-terms as pertur- 
bation. 

When J' = 0, all type-A tetrahedra are decoupled 
from each other, and we can easily construct eigenfunc- 
tions of the entire system from those of one tetrahedron. 
On an isolated tetrahedron, eigenfunctions of its total 
spin Sa are also eigenfunctions of the Hamiltonian. Since 
there are four spins each with S = ^, Sa is either 0, 1, or 
2, with degeneracy 2, 3, or 1, and energy — | J, —^J, or 
+ ^J, respectively. In the ground-state multiplet, each 
type-A tetrahedron takes one of the two 5a=0 states. 
Therefore, the degeneracy is 2^* {Na- the number of 
type-A tetrahedra, i.e., a quarter of the number of sites), 
and all the states are of course spin singlet. The lowest 
excited multiplet is a spin triplet, and the spin gap is J. 

As we will see later, it is convenient to use a chiral 
basis for each tetrahedron. In the S^ = subspace, two 
spin-singlet states are 

|0, +) = [(12)+ + e^*(13)+ + e'»ni4)+]/V3, (2) 
|0, -) = [(12)+ + e-**(13)+ + e'^(14)+]/V3, 

where (*j)± = (^-^'l TTTT) ± S+S+\ iiii))/V2 and 
(/)=27r/3. As for the spin triplet, we obtain two chiral 
states, |1,+), |1,— ), by replacing with {ij)-, com- 

pleted by the non-chiral state |1,0) = [(12)_ + (13)_ + 
(14)_]/^/3. The spin-quintet state is also non-chiral, 
|2,0) = [(12)+ + (13)+ + (14)+]/V3. Let us define the 
spin chirality by Xijk — "^^i ' i^j ^ ^k), and these bases 
are also the eigenfunctions of the chirality around site 1, 
X234|'S', r) = t\S,t) (r — +,—,0). We can also define 
the chirality around other sites in a similar way, and it 
is important that the spin singlets do not change their 
chirality, while the triplets do. This is easily understood 
since the singlets belong to the complex ID representa- 
tions, r2 and Fa, of the tetrahedral point group, whereas 
the triplets belong to the 3D representation T2. 
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Fig. 3. Supertetrahedron comprised of four type-A tetrahedra, 
and value of parameter Oj in the effective model. 



When the J'-term is switched on, this leads to a mix- 
ing of different states in the ground-state multiplet, and 
we may expect a lift of the thermodynamic degeneracy. 
First, consider a pair of neighboring type-A tetrahedra, 
and connect them by a single bond with exchange cou- 
pling J'. I show in Fig. 2(a) the energy of 70 states in 
S'^Qt = subspace, as a function of J'. The degeneracy 
of the ground-state multiplet is 4-fold, and it is never 
lifted. Therefore, we have to go to the next step, and 
consider three neighboring type-A tetrahedra each pair 
of which is connected by one bond. There are two differ- 
ent cases. In the first case, two perturbative bonds are 
connected to the same site on each tetrahedron, whereas 
in the second case they are connected to two distinct 
sites. It is found that the second case lifts the degener- 
acy, but the first case does not. Figure 2(b) shows the 
energy of 924 states in the S^^^ = subspace, for the 
second case. The ground-state multiplet with 8-fold de- 
generacy is now split into 4 levels with 1, 3, 1, and 3 fold 
degeneracies, respectively, in the ascending order of en- 
ergy. At J' = J, this energy split is 0.2468 J between the 
lowest and highest levels, while the singlet-triplet gap is 
0.5491J. 

These behaviors are described by degenerate pertur- 
bation in J'. The second order only yields a constant 
self-energy, — ^J'^/J per spin, and we require third- 
order terms to mix different chirality configurations on 
a tetrahedron triple. To describe chirality for the tetra- 
hedron positioned at R, I introduce spin-^ pseudospin 
operators, ^tj^, such that their z-component is a diago- 
nal matrix in the chirality space. Using these operators, 
the third-order effective model is written aa3' 

H^l^^ E [c-&(l-2e(ai).rR) (3) 
(R,R',R"> 

x(l - 2e(a2) • rR,)(l - M^^s) ' ^R^O], 

where c= ^.J'^ / J^, b = J'^/ J^, e(Q;) = (cos a,sin a,0) 
Type-A tetrahedra constitute the fee arrangement, and 
the summation is taken over all the tetrahedron triples 
of the second case explained above. Note that the to- 
tal chirality is not conserved and also that the chirality 
pseudospins interact only in their transverse directions. 
The angle aj depends on the position of not only its cor- 
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Fig. 4. Example of virtual processes in (a) the second and (b) the 
third order in J'/ J. Balls are type-A tetrahedra with chirality 
(r=+, — ,0) indicated around each bond. Perturbation in each 
step is indicated by a solid bond. 



responding tetrahedron but also the other two in each 
triple. Its value is shown in Fig. 3. The one-body terms 
in eq. (||) are canceled out after taking the summation. 

We may understand the second and third-order effec- 
tive models by examining their virtual processes. It is 
essential to recall that the chirality of each spin-singlet 
eigenfunction on a tetrahedron is identical around any 
of its four sites. On the other hand, a spin-triplet chiral 
eigenstate around one site is a mixture of all the three 
chiral eigenstates around another site. It is also impor- 
tant that perturbation of a single bond conserves the 
chirality of each tetrahedron around the site to which 
the bond is connected, in both spin singlet and triplet 
multiplets. Therefore, as shown in Fig. 4(a), the chiral- 
ity configuration on a tetrahedron pair always returns to 
the original one, after the second-order virtual processes, 
more generally in any orders, and the ground state de- 
generacy is not lifted. In the case of a coupled tetra- 
hedron triple shown in Fig. 4(b), chirality is mixed up 
in the virtual processes including spin-triplet tetrahedra, 
since two bonds are connected to different sites on each 
tetrahedron. Therefore, any of the eight chirality config- 
urations can undergo a transition to all of them. This 
is why the effective model (||) has one-, two-, and three- 
body interactions in the transverse components in the 
chirality space. 

Now I discuss the ground state of the effective model 
(^) using the mean-field approach. To this end, I first 
determine the mean-field ground state of its elementary 
building block shown in Fig. 3, the "supertetrahedron" , 
each of its vertices (j=l'^4) being a type-A tetrahedron. 
The mean-field Hamiltonian, E, is simply obtained by 
replacing chirality operators in eq. (||) with classical unit 
vectors, t (cos0, sinfl, 0), and aside from a constant 
it is given as 

£{61,92,93,64) = -6{(1 - 2 cos 02) [(1 - 2cos(03 + 0)] 

X [(1 - 2 COS(04 - 0)] + {02,03, 94^91,64,93) 

+ (^2, 03-S4 ^ 04, 01,02) + {02, 03,04 03, 02, 0l)} ■ 

(4) 

This has a symmetry associated with spatial rotation: 

£{64,62,63,04)^ E{01 + ^,03 + <f>, 04 + ^,02 + (!>)= E{94 + 
,P,02 + ^,0l + <j>, 03 + (iy)=E{92 + <l),94 + <l^, 03 + <^,^l + </-) = 

E{03+<t>, 0i+<t>, 02+<t>, 04+(f)- I minimize this classical en- 




Fig. 5. Mean-field solution of the effective model in the chirality 
space. Azimuthal angle 9 of the ordered pseudospin is shown. 



ergy and find four types of ground-state configurations, 

(01,^2,03,^4)= (0i,^,i^,|7r), (^,02,§7r,i^), (5) 
(iyr, |7r,03,7r), (|7r, ^tt, tt, ^4). 

In each solution, one angle variable is still free and the 
classical energy is degenerate with respect to it. 

I can patch up these four types of solutions for a su- 
pertetrahedron to construct a ground-state configuration 
on the entire fee lattice, and the result is shown in Fig. 5. 
Three of the four sublattices show a long-range order. 
The other one remains disordered, and its residual en- 
tropy is log 2 per original spin. I emphasize here that 
this solution is the global energy minimum within the 
mean-field approximation. This is because the solution 
minimizes the energy for each supertetrahedron in the 
whole fee lattice without any extra energy cost. The de- 
generacy of this state is 4 x 2^*/^, where the first factor 
is from the number of sublattices and the second one is 
from the disordered sublattice. 

The long-range order contributes to the chirality struc- 
ture factor a set of Bragg peaks at the reciprocal lattice 
vectors of the simple cubic lattice, Q = 2TT{nx,ny,nz) 
(n's are integer), in addition to a constant background 
due to the disordered sublattice. The amplitude of the 
Bragg peaks is given by S'^(Q)=A^a^Er.R' e*Q (^-^') 
('^R'^r) =|(Caz+C2a;+C:ry), wheie Cy^ = sin^(i((3y- 
Qz))— |(1 — (—1)"!'+"^) and so on, and therefore van- 
ishes at Q where n^, riy, and have the same parity. 

The mean-field ground state has a finite expectation 
value of the pseudospin, ("rp^) — (cos0, sin0, 0). How- 
ever, its z-component describing spin chirality is zero, 
meaning that the time-reversal symmetry is not broken 
and dominant spin-spin correlations are coUinear rather 
than chiral. The transverse components are finite in- 
stead, and I now discuss its physics with considering the 
wavefunction of four spins on an original tetrahedron. 
This for the given value of angle 6 is given as the eigen- 
function of the operator, (x) ■ t — cos 9 + sin 9. 
All of the three cases, 6—n,^Tr,^TT, correspond to a va- 
lence bond state, i.e., a dimer pair, but with different 
pair combinations: 

= 7r : [12](»[34], ^tt : [13]®[24], fjr : [14]®[23], (6) 
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where [m'] denotes the spin singlet wavefunction of two 
spins, i and i' (see Fig. 1). Therefore, the mean-field 
ground state has a long-range order of dimerization in 
three of the four tetrahedron sublattices and breaks the 
translation, mirror, and rotation symmetries in space. 
However, it does not break symmetry in spin space, since 
the order occurs within the spin-singlet multiplet. The 
possibility of dimerization was first discussed by Harris et 
al.0 and the present solution agrees with their proposal. 

Koga and Kawakamia^ also discussed dimerization on 
the pyrochlore lattice with including additional distor- 
tion which favors the q = order of either a tetramer 
singlet or one of the dimer pairs in eq. (^, and found 
that these two states become degenerate as the distor- 
tion is switched off. Their solution is not continuously 
connected to the present solution, but this is because 
the distortion under consideration is limited to intra- 
tetrahedron modes with q = 0. As the result of Sr{Q) 
shows, the present solution has a triple-q structure with 
wavevectors (27r,0,0), (0,27r,0), and (0,0,27r). The dimer- 
ized state resonates with the tetramer state to realize 
the other two dimerized states in eq. (^) with spatially 
modulated coefficients, and the three-body interactions 
in the effective model (||) stabilize the triple-q structure 
with four sublattices. 

To compare with the experiments, one useful quan- 
tity is the spin structure factor, S'o-(q)= {4:NA)~^J2r v 
g»q (r-r') . g^^)^ rpj^j^ calculated for the mean-field 
ground state and the result is 

Sa{(l)^j[3~ (cos cos cos ig^ cos ig^, (7) 
+ cosjqxCOSjqy)]. 

Since this state is a spin singlet with no magnetic long 
range order, the spin structure factor has no Bragg peak 
and continuously changes with the wavevector. It has a 
maximum value at q's where one of cos {iJ,=x,y,z) 
is and another one is —1, e.g., Qi = (47r, 47r, 0) and 
Q2=(47r, 27r, 0). Therefore, these q's form three orthog- 
onal sets of parallel lines, reflecting strong geometrical 
frustration. This result (|^) is qualitatively similar to the 
result. Fig. 4, of Canals and LacroixjZP but the maximum 
position is shifted as for Qi in their calculation. On the 
one hand, this may be because spins on different tetra- 
hedra have no correlations within the present approach. 
On the other hand. Canals and Lacroix used the second- 
order perturbation in the inter-tetrahedron interaction 
J'. As I pointed out above, up to the second order, the 
J'-term induces no spin-spin correlations between dif- 
ferent tetrahedra for the singlet ground-state multiplet. 
Therefore, the finite shift from Qi in their result is due 
to the contribution of excited-state multiplets at finite 
temperatures, and it will vanish with the temperature 
approaching zero. We could improve the present result 
by also including perturbation for the wavefunction, but 
it is beyond the scope of this paper. 

In this paper, I have examined the ground state of 
the 5=^ Heisenberg antiferroniagnet on the pyrochlore 
lattice. Assuming the ground state is a spin liquid, I 
obtained the effective model for a thermodynamic num- 
ber of spin-singlet states below the spin gap and exam- 



ined symmetry breaking, based on the third-order de- 
generate perturbation starting from decoupled tetrahe- 
dra. I determined its ground state using the mean-field 
approximation. The result is a long-range order of spin- 
singlet dimers with a triple-q structure, and three of the 
four constituent tetrahedron sublattices have a ordered 
dimer-pair pattern. Therefore, spin-triplet excitations 
have a finite energy gap, whereas spin-singlet excitations 
exist within the spin gap. These singlet excitations are 
oscillations of the dimerization pattern in the ordered 
sublattices. One tetrahedron sublattice remains disor- 
dered, and it is randomly dimerized, tetramerized, even 
chiral, or any superposition of these. We expect that this 
sublattice will ultimately show a long-range order once 
quantum fluctuations are taken into account beyond the 
present approach. In a separate paper, I will discuss this 
problem both for spin and chirality degrees of freedom. 

I have investigated in the present study the lift of 
thermodynamic degeneracy driven by spin exchange pro- 
cesses, but there may exist oither types of instability. Re- 
cently, Yamashita and Uedallj) found that a Jahn- Teller 
lattice distortion lifts this degeneracy for the 5=1 anti- 
ferromagnet, and they proposed that this is a mechanipri 
of a structural phase transition observed in ZnV2 04.M' 
This compound also undergoes a magnetic phase tran- 
sition at a lower temperature. To explain its magnetic 
order, theoretical approaches limited to the spin singlet 
multiplet are not sufficient and we must include finite- 
spin multiplets not only in the virtual processes but also 
in the real processes. This problem is left for a future 
project. 

I thank Frederic Mila for valuable discussions. This 
work was supported by a Grant-in- Aid from the Ministry 
of Education, Science, Sports and Culture and also by 
University of Tsukuba Research Projects. 
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